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Abstract 

Let v be a vertex of a graph G. By the local complementation of G at v 
we mean to complement the subgraph induced by the neighbors of v. This 
operator can be generalized as follows. Assume that, each edge of G has a label 
in the finite field F g . Let (gij) be set of labels (gij is the label of edge ij). 
We define two types of operators. For the first one, let v be a vertex of G and 
a € F q , and obtain the graph with labels g'^ = g^ + ag„ig V j. For the second, 
if 7^ b G F q the resulted graph is a graph with labels g'^ = bg V i and g'^ = g^ , 
for i,j unequal to v. It is clear that if the field is binary, the operators are just 
local complementations that we described. 

The problem of whether two graphs are equivalent under local complemen- 
tations has been studied, [_>]. Here we consider the general case and assuming 
that q is odd, present the first known efficient algorithm to verify whether two 
graphs arc locally equivalent or not. 

1 Introduction 

A labeled graph is a graph all of whose edges have a label chosen from a (finite) field. 
This definition covers the usual graphs when one restricts the filed to the binary 
field, F2. We want to define the notion of local equivalency over (labeled) graphs 
but, for simplicity, let us first consider the binary case. In the binary case, i.e., 
when the field is F2, consider the following operation, called local complementation. 
Choose a vertex, and replace the subgraph induced on the neighbors of this vertex 
by its complement. Two graphs are called locally equivalent if one can obtain one 
of them from the other by applying some local operations as above. 

In general, when the field is not binary, two types of operators are involved. The 
first one is just the generalized version of the operator in the binary case. Let the 
graph G be labeled with labels forming a symmetric matrix G = {gij) with zero 
diagonal over F q , where q is a power of an odd prime number, and F q is the field 
with q elements. Let v be a vertex of this graph, and a € F q . We define G * a v 



'Massachusetts Institute of Technology, Mathematics department & Computer Science and Ar- 
tificial Intelligence Laboratory, m_bahram@mit.edu. 

^Massachusetts Institute of Technology, Mathematics department, salman@math.mit.edu. 



1 



to be the graph with labels G' = (g[j), where g'^ = g, t j + ag v ig v j. In the second 
type of operators we multiply the edges incident to a vertex v by a non-zero b € F 9 , 
and denote this graph by G o b v. In other words, G o b v is the graph with labels 
G" = (g'ij), where g" { = bg v i and </y = for i,j unequal to v. Similar to the 
previous situation, two graphs are called locally equivalent if one of them can be 
obtained from the other by applying a series of operators * and o. 

Studying and investigating the local equivalency of graphs has become a natural 
problem in quantum computing, and playing a significant role especially in quantum 
error correcting codes, due to the recent work of [I], [2], [6] and [7] . Namely, in the 
quantum computing setting, some states, called graph states, have a description as 
the common eigenvectors of a subgroup of the Pauli group. These states are called 
graphs states because their associate subgroup is defined bases on a labeled graphs. 
Using graph states, we may be able to create more preferable quantum codes, due 
to the property that the obtained codes have relatively shorter descriptions, and 
are more algebraically structured. Hence, combining the theory of quantum error 
correcting codes and the tools in graph theory, leads us to describe and investigate 
the properties of graph states more and more deeply. 

The key point is that, we can obtain one graph state from another by applying 
elements of, what is called local Clifford group. If two states are equivalent under 
local Clifford group, they present similar properties in quantum computing. In fact, 
as shown in [2] and [6] , two graph states are equivalent under the local Clifford group 
if their associated graphs are locally equivalent by the local operators described 
earlier. So, this question is coming up naturally that, when two graphs are equivalent 
up to these operators, and how we can recognize them. 

The special case of q = 2, has been studied in the work of Andre Bouchet, [3J, 
[5], and a polynomial time algorithm for recognizing the local equivalency of two 
(simple) graphs is described in [3j. In fact, he showed that, for any two graphs there 
is a system of equations such that, the two graphs are locally equivalent iff those 
equations have a solution. 

Same as binary case, when q is odd, recognizing the locally equivalent graphs is 
equivalent to solving a system of equations, some of which are linear and the rest 
are quadratic. But, their algebraic structure are different, and is more non-linear 
compared to the binary case. Indeed, in the binary field every element satisfies 
a 2 = a, and hence quadratic equations on binary field exhibit linear properties. 
The algorithm described in [3] takes the advantage of this property of F2. The 
situation in non-binary case is completely different, and the quadratic equations do 
not exhibit linear properties in general. In the present paper, we study in details 
the structure of the solutions of these equations, and present an efficient algorithm 
to solve the problem of recognizing locally equivalent graphs. 

1.1 Main ideas 

The main ideas in this paper are as follows. First of all, we introduce isotropic 
systems that are geometrically known objects, and define an equivalency relation on 
them. In this definition two isotropic systems are equivalent if a system of equations 
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has a solution. Then, we define the isotropic system associated to a (labeled) graph, 
and show that two graphs are locally equivalent if their associated isotropic systems 
are locally equivalent. Using this idea, we convert the problem of local equivalency 
of graphs to the existence of a solution for a system of algebraic equations. 

Unfortunately, these equations are not all linear. So in general, it is hard to 
decide whether there is a solution or not. But, in our case, their solutions have 
some nice properties. In fact, we prove that, if there is one solution then there are 
many. In other words, if two graphs are locally equivalent then, in some sense, there 
are many solutions for their associated system of equations. 

The idea to prove this property, is to correspond the solutions of the system 
of equations for two graphs, to the solutions to some equations associated to just 
one of them. What we call them internal solutions. In fact, we show that, instead 
of studying local operators that convert one graph to the other, it is sufficient to 
know set of local operators that send a graph back to itself. This correspondence 
allows us to somehow consider a linear structure for the set of solutions, and to show 
that if there is a solution then the solutions contain an affine subspace of constant 
co-dimension. In other words, if there is a solution then there are many, and so it 
is not hard to find one of them. 

2 Isotropic systems and locally equivalent graphs 

Assume that p is an odd prime number, and F g is the field of q elements with 
characteristic p. Suppose that G is a graph on n vertices. We call G a labeled graph 
on F q , if labels of its edges form an n x n symmetric matrix G = (gij) with zero 
diagonal over F q , where is the label of edge ij. 

2.1 Local operators over graphs 

Definition 2.1 Let G be a labeled graph with labels forming a symmetric matrix 
G = (gij) on F q . For vertex v of G and a £ F q , define G * a v to be a graph with the 
label matrix G' = (g^) such that for all i, g' vi = g v i, and for i,j unequal to v, 

9ij = 9ij +ag vi g vj , 

and moreover, G' is symmetric with zero diagonal. 

Also, for a non-zero number b € F g define G o b v to be a graph with the label 
matrix G' = (g^) such that for all i, g' iv = bgi v , and g[ - = g^ if i,j are unequal to 
v, and again, G' is symmetric with zero diagonal. 

Two graphs G and G' are called locally equivalent if there exists a sequence 
of above operations that acting on G gives G' . Notice that, these operations are 
invertible, so that this is an equivalency relation. 



3 



2.2 Isotropic systems 

Let Fg be the n-dimensional vector space over F q , and consider the standard bilinear 
form on it. That is, for vectors 

X = (xi,-- - ,x n ), Y = (yi,--- ,y„), 

in F™, define 

n 

(X,F)=^x i y i . 
i=l 

( ., .) is a non-degenerate, symmetric bilinear form. Using this form, we define a 
non-degenerate anti-symmetric bilinear form on V = F^™, the 2n-dimensional vector 
space over F q . For vectors (X, X') and (Y, Y') in V, set 

((X,X'),(Y,Y')) = (X,X').A-(Y,Y') T , 

where 




In other words, 

((X,X%(Y,Y')) = (X,Y')-(X',Y). 

Due to the nature of anti-symmetric forms, we are now in the situation of introducing 
a geometrically known concept, isotropic systems. 

Definition 2.2 A subspace W of V is called an isotropic system , if it is an n- 

dimensional subspace and ((X, X'), (Y, Y')) = 0, for any (X, X'), (Y,Y r ) G W. In 
fact, since dimW = n and (.,.) is a non- degenerate bilinear form, we have 

W = {V G V : (V,W) = 0, VW G W}. (1) 

In this paper the basic examples of isotropic systems are isotropic systems as- 
sociated to graphs. For every graph G, let Wg be the vector subspace generated by 
the rows of matrix (/ j G) . That is 

W G = {X ■ (I | G):le F£}, 

where (/ | G) is a matrix with two n x n blocks which the first one is identity. 
Since, G is a symmetric matrix 

(/ | G)-A-(G | -I) T = G T -G = 0, 

and we conclude that Wg forms an isotropic system which is called the isotropic 
system associated to G. 



2.3 Locally equivalent isotropic systems 



Definition 2.3 Suppose that A is a 2n x 2n matrix 



A 



( 



Z T 
X Y 



consisting of four diagonal matrices X = diag(x\, - ■ ■ ,x n ),Y = diag(y\, • • • ,y n ), 
Z = diag(z±, • • • , z n ) and T = diag(t±, ■ ■ ■ ,t n ). A is called normal if 



Notice that, every normal matrix is invertible, its inverse is normal as well, and in 
addition the multiplication of normal matrices is again normal. In particular, when 
Z = I and X = 0, we call A trivial. 

For an isotropic system W and a normal matrix A, define 



Since A is normal, one can verify that VV • A is also an isotropic system. 

Two isotropic systems W and W' are called locally equivalent if there exists a 
normal matrix A such that W" = W • A. By the properties of normal matrices, it is 
clear that this is an equivalency relation. The importance of this relation would be 
clear once we state the following theorem. 

Theorem 2.1 Two graphs G and H on the same vertex sets are locally equivalent 
if and only if their associated isotropic systems are locally equivalent. 

Proof: For the only if part, it is sufficient to show that the systems WcWgv, 
and also the systems Wg, Wgojj are locally equivalent. First, note that the rows of 
(/ | G * a i) form a basis for Wc* a i- Let 



where X = diag(0, • • • ,0, —a, 0, • • • , 0), Y = I, Z = I and T = 0. We just need to 



j ^ i, and k ^ i, the product of these two rows is (gjk + cigikgij — agijgik) — (djk) = 0. 
If j = i and k ^ i, it is (g-ik) — (g-ik) = 0, and finally, if j / i, and k = i then this 
product is again (gij) — (gij) = 0. Also, the i-th row of these matrices are equal, and 
therefore, according to the definition of the inner product, these rows are orthogonal. 
Thus, Wg and Wc* a i are locally equivalent. 
Similarly, for Wc,WGo b i, let 



Y ■ Z — X T = I. 



W • A = {W ■ A : W £ W}. 
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where X' = 0, T' = and 



Y' = diag(l, ■ ■ ■ ,1,6,1,- ■■ , 1), Z' = diag(l, ■ ■ ■ ,1). 



Similar to the previous case, one can easily check that the rows of (/ | G) • B and 
(/ G o b i) are orthogonal, and therefore, Wg, yVco b i are locally equivalent. 

To prove the if part, suppose that Wg* and Wh are locally equivalent. Then 
there exists a normal matrix 

V. 

where, 



X = diag(x!,--- ,x n ),Y = diag(y lr -- ,y n ), 

Z = diag(z!,- ■ ■ ,z n ),T = diag(h,- ■ ■ ,t n ), 

and rows of (/ | G) A form a basis for the isotropic system Wh- Therefore, there 
exists an invertible matrix U such that U ■ (/ i G) ■ A = (/ | H^j . 
For every i, let 

'Zi T- 

. Xi Yi 



Ai 



where 



l\Ti 



diag(0, 



AU,o, 



a) 



Zi = diag(l,--- ,1,^,1,> 

Xi = diag(0, ■ ■ ■ , 0, Xi, 0, • • ■ , 0), Yi = diag(l, ■ 

The matrices ^4j's all commute and A = A\A<i - ■ ■ A n . 

We prove the theorem by induction on the number of non-trivial matrices A^. If 
all AiS are trivial then A is also trivial, and we have 

(/ | G)-A=(I | G)-(J ^) = (/ | T + GY). 

Therefore, we have U ■ (I \ T + GY) = (/ | -ff) . Looking at the first blocks in 
this equation, we get that U = I, and T + GY = H. The diagonal of H is zero, and 
hence T = 0. Also, A is a normal matrix, so that, A = I and G = H. 

Hence, suppose that at least one of A^s is non-trivial. We consider two cases: 

Case(i). Zi / for some zq, where Ai is non-trivial. Let (J j G) Ai = 
(V | D). Then, 



/I 
1 



V = 



\0 



o \ 





'•■ 

■■ 1 J 
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In order to get to the inverse of V, we should multiply the io-th row by z^ 1 , and then 
by —Xi gji and add it to the j-th row, for any j ^ iq. Thus, V^ 1 (/ ; G) Ai = 
(/ i y _1 Z)), and the jk-th entry of V~ 1 D, for j unequal to k, and io, 



is 



(V- l D) ]k 



9jk 



Z i x io9ioj9i k- 



Also, for j ^ i , 



9k) j ■ 



The matrix V~ l D may have non-zero entries on its diagonal. But, by elemen- 
tary linear algebra, there exists a trivial matrix A', such that all the entries of 
V^ 1 (/ | G) Ai A' are equal to V^ 1 (I j G) A{ except those, on the diagonal 
of the second block, which are zero. Therefore, the rows of V^ 1 (/ j G) A^A' 
span an isotropic system associated to some graph, and by the above equalities, this 
graph is nothing but G *, -i s i o, -u i . 

z i ( z i > 

On the other hand, we have 

(/ | H)=U(I | G) A = UV{V- 1 (I | G) A^A'^A'^A"), 

where, A" is equal to the multiplication of all Aj's, except A io . 

Now, V^ 1 (/ \ G) Ai A' is an isotropic system, associated to the graph 

G* { _ z -i Xio) io ° ZlQ io- 

Also, the number of non-trivial terms in A'~ l A" is strictly less than the number of 
non-trivial terms in A, and therefore, by induction, we obtain the desired result. 

Case(ii). Zi = for all i's, where Ai is non-trivial. Notice that, in this case xi ^ 
for any i, where Ai is non-trivial. Because, A is a normal matrix and yiZi — Xit% — 1. 
Suppose that, Ai is non-trivial. If for every non-trivial Aj, gi j = 0, then the io-th 
row of the first block of (/ j G) A is zero. Hence, it is not invertible and the first 
block of U (I j G) A can not be identity. Thus, there exists an ii, such that A^ 
is non-trivial and g^ ^ 0. Therefore, the first block of (l j G) A^A^ is 



V = 



(\ 









x io9lio •EiiSlii 



o \ 




x ii9ioii 



•Eio9nio x i\9ni\ 



J 



In order to invert V, one has to multiply the io-th and the ii-st rows by 
{ x h9ioiiY l an d (x^g^)' 1 respectively, and then multiply the row ii by -x io g jio 
and add it to the j-th row, for any j. Also, exactly the same process for row io must 
be done. After this, we should change the rows io and i\. By this process, we get 
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to a matrix with identity in the first block and a symmetric matrix on the second 
block. But, the diagonal entries of this block may be non-zero. In order to handle 
this issue, by multiplying by an appropriate trivial matrix A', we get 

V- 1 (/ | G) A i0 A h A' = (I | G') , 

where G' = (g^) is the graph with entries 

9jk = 9jk — 9i$i x 9ioi9i\k — 9i i 1 9iij9iok, 



9i j — x i Qioi^iiji 



Jioh x i x i 1 5i ii" 



for all j, k ^ iq, ii, and 

for all j ^ io, ij. Moreover 
Therefore, 

G ' = G i0 H i0 h H i0 l ° H - 

We have (I \ H) = UV (I \ G') A'- 1 A", where A" is equal to the multi- 
plication of all Aj's, except Ai and Ai x . Also the number of non-trivial terms in 
A'^A" is strictly less than this number in A, and therefore, by induction, the result 
is proved. 

□ 



3 System of equations and normal matrices 

Let G and H be two graphs, and g, h be their neighborhood functions respectively, 
meaning that g(i) is a vector such that its j-th coordinate, gij, is the label of the 
edge ij in G. The same holds for the graph H and the function h. Using theorem 
12.11 and relation (pQ) in the definition of isotropic systems, one obtains that G and 
H are locally equivalent if there exists a normal matrix 

A= (x y) ' 

such that rows of (/ j G) ■ A are all orthogonal to the rows of (/ | H^j . This 
condition is equivalent to the following: 

(X,g(i) x h(j)) - (Y,g(i) x ej ) + (Z,* x h(j)) - (T,e, x ej ) = 0, (2) 

for any two vertices Here, by u x v for two vectors u and v of the same size, 
we mean a vector of the same size whose A;-th coordinate is the product of the k-th 
coordinates of u and v. Also, is a vector whose all coordinates are zero except 
the i-th, which is one. Notice that, in this formula, and also later on, we look at the 
diagonal n x n matrices as vectors of size n when necessary. 
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Let F^ n be the space of the vectors of the form (X, Y, Z,T), provided with the 
following symmetric bilinear form 

((X, Y, Z, T), (X',Y', Z', T')) = (X, X') - (Y, Y') + (Z, Z') - (T, T'). 

Moreover, for any pair of vertices we define the following function, so called 
Lambda Function, which plays a key role in the whole section. 

Hhj) = \9(i) x h(j),g(i) x ej,e» x h(J),ei x eA 

and 

X(G,H) = Span{X(i,j) : 

Notice that, equation ([2]) is equivalent to say that X(i, j) is orthogonal to (X, Y, Z, T). 
Therefore, the problem of local equivalency of G and H reduces to the following: 

Two graphs G and H are locally equivalent iff there exists a vector (X, Y, Z, T) 
orthogonal to X(G, H), and Y X Z — X X T = I. 

3.1 Changes of X(G,H) under local operators 

We develop some more technical tools to attack this problem. For <p = (X, Y, Z, T) 
G Ff , let 

^ = (-Z,-T,-X,-Y), 

^ 2 = (y,x,r, z), 

and for a G F^ n , a G F q and I = 1,2 define 

<fi *ia a = (j) — a<p l x a. 

Also, for any subspace iV of F^ n , let 

N *i, a a = {4> *i A a: 4>e N}. 

Lemma 3.1 Let a G F^ n be such that a x a 1 = for I = 1,2. For <f>, ip G F^ n , and 
any subspace N ofFg n , the following properties hold: 

(i) 4> — > 4> *i,a a is bijective. 

(fi) (4>*l,aU,1p *l,aU l ) = {(j),i>). 

(iii) (Nn, a a) L = N L n,aa l . 

Proof: For (i), by a simple induction one can check that after k times iteration, 
we get 4> *i a a ■ ■ ■ a = <f> — ka<p l x a. Therefore, for k = p we end up with the 
identity map, and hence <p — > <p *i,a a is a bijection. 
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For (ii), using the facts that (</> x ip' l ,ijj l ) = —(<P l x V'' an d ((f>,i> x i>') = 
(0 x ip',ip), we have 

(0 *i :a a, if) *i >a a 1 ) = (4>, ip) — a{4>, ip l x a 1 ) — a((f) 1 x a, ip) + a 2 ((f) 1 xa^'xa 1 ) 
= (</>,-!/>) - a(4> x a', - x a,ip) + a 2 (4> 1 x a x a',-*//) 

(iii) is a direct consequence of (ii). 

□ 

In the next two theorems, we study the effect of local operations on the set 
X(G,H)- L , and observe that this set is well-behaved under these types of operators. 

Theorem 3.1 Let G and H be two graphs on the same vertex set with neighborhood 
functions g and h, respectively. For every vertex i, 

(i) X(G * a i, H) = X(G, H) *i, a {-g(i) x g(i), -g(i) x g(i), - ei , -a). 

(ii) X(G, H * a i) = X(G, H) * 2 , a (h(i) x h(i), a, h(i) x h(i),&i). 
And more importantly, 

(iii) X(G * a i, H) L = X(G, H) x * t , a (e;, g(i) x g(i),g(i) x g(i)). 

(iv) X(G, H * a i) 1 = X(G, H) 1 * 2 , a (ei, fc(t) x h(i),e u h(i) x h{i)). 

Proof: The proof of (ii) is similar to that of (i). Also, parts (iii) and (iv) are 
immediate consequences of (i), (ii) and the third part of lemma I3TT1 Hence, we just 
need to prove the first part. 

Suppose that g' is the neighborhood function of G * a i. Then, we have 

g'{j) = 9(j) + a 9ijg{i) - agfjej. 
The space X(G * a i, H) is generated by A'(j, k)'s, where 



X'(j,k) = (g'(J) x h(k),g'(j) x e k ,ej x h(k),ej x e k 

= X(j,k) +agij(g(i) x h(k),g(i) x e fc ,0,o) -ag^ej x h(k),ej x e fc ,0,0 
= X(j, k) + agij (a(i, k) - (0, 0, a x h(k),ei x e fc )) 
-agfjUj x h(k), ej x e fc ,0,o). 

Set X"(j, k) = X'(j, k)—agijX'(i, k). It is easy to check that X"(J, k)'s also generate 
the space X(G * a i, H), and we have 
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□ 



A"(j, k) = A(j, k) + agij\(i, k) - agij (o, 0, e { x h(k), x e fe ) 
-agfj^ej x h(k),ej x e fc ,0,o) -agijX(i,k) 

= X(j,k) -a(o,0,g(j) x a x h(k),g(j) x e { x e fc ) 

-a(g{i) x g(i) x ej x h(k),g(i) x g(i) x e,- x e fc ,0,o) 

= A(j, fc) -a(ej x h(k),ej x e fc ,g(j) x h(k),g(J) x e fc ) 

x (5(0 x 9(i),9(i) x ej, 
= A(j, fc) *i j0 ( - fli(i) x g(i), -g(i) x g(i), - ei , -e») . 

Therefore, A(G * a i, #) = A(G, ff) * 1>a (-g(i) x x #(T), -e,, -e,). 

To state the next theorem define 

h,i = I +{b-l)ei = {!,■■■ ,1,6,1,- •■ ,1), 
for every 7^ 6 G F g and vertex i. 

Theorem 3.2 Let G and H be two graphs on the same vertex sets, and with neigh- 
borhood functions g and h, respectively. For every vertex i, 

(i) The map <j) -> <j) x (f bui , f b2>i , f b:iti , f bi>i ) is bijective, for non-zero elements 
bi,b 2 ,b 3 ,b 4: £ F q . 

(ii) X(Go b i,H) ± = X(G,H) ± x (f b -i ti ,f b -i,i J b ,i ,f b ,i). 

(iii) A(G, H o b i) 1 - = A(G, H)^ x (f b . lti , / M , f b . lti , 

Proof: The proof of (i) is straight forward. To prove (ii), notice that g' , the 
neighborhood function of G o b i, is given by g'(j) = g(j) x f b) i if j 7^ i, and g'(i) = 
bg(i). Hence, if (4>,X(j,k)} = 0, then 

{<P x (h-ij ,f b -i,i ,f b ,i ,fb,i),\'(j,k)) = 0, 
where \'(j,k) = (g'{j) x h(k),g'(j) x 

e k 1 e j X h(k), ej x e^j . Part (iii) is similar to 

(ii). 

□ 

3.2 Changes of determinant function 

We now define the determinant function for a vector in F 4n . For a vector 4> = 
(X,Y,Z,T), set 

det(f> = YxZ -XxT. 
Some straight forward computations easily lead to the proof of the following lemma. 
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Lemma 3.2 



(i) detcj) 



detcj)* 1>a ((%,(%, g(i) x g(i),g(i) x g(i)). 



(ii) detcj) 



det cf> * 



2,a (ei,g(i) x g(i),ei,g(i) x g(i)). 



(iii) det 



det<?!> x (fb-i t i,f b -i ti Jb,i,fb,i)- 
detcj) x (f b -i,i,fb,i,fb-i,iifb,i)- 



(iv) det 



□ 



Note that, the functions that we see on the right hand side of parts (i) to (iv) in 
this lemma, are exactly the ones appear in theorems 13. II and 13. 2[ Hence, this lemma 
states that the determinant function is invariant under the action of * and o. 

3.3 What is the new picture? 

In the new setting, the problem of verifying whether or not two graphs G, H are 
locally equivalent, is equivalent to finding a vector cp E F^ n such that cf> G X(G,H)- L 
and det cj> = I. 

To get a more convenient notation, let A(G,H) = X(G,H) ± , and a(G,H) be 
the set of solutions, i.e., vectors cp £ A{G,H) satisfying detcj) = I. Then, we get to 
the following picture from theorems I3.1|. 13.21 together with lemma 13.21 

If graphs G\,G 2 are locally equivalent, as well as the graphs Hi,H2, then there 
exists a (linear) bijection (3, such that 



Even though the function (3 depends on Gi,G2,H\ and H2, but it gives us useful 
information on the locally equivalent graphs. Namely, if two graphs G and H are 
locally equivalent, then roughly speaking, the relative linear position of cr(G, H) 
inside A(G, H) is exactly the same as the relative linear position of cr(G, G) inside 
A(G, G). For instance, once we prove that for every graph G, a(G, G) is a large sub- 
set of A(G, G), in the sense that it contains a linear subspace of small co-dimension, 
then the same must be true for a(G, H) inside A(G, H), when G and H are locally 
equivalent. 

4 Internal solutions 

We have shown that, for locally equivalent graphs G\ and G2, and again locally 
equivalent graphs H\ and H2, there exists a linear bijection /3, such that 



A(G 2 ,tf 2 ) 



P{k{Gx,H x )), 



a(G 2 ,H 2 ) 



f3(a(G 1 ,H 1 )). 



A(G 2 ,H 2 )=P(A(G 1 ,H 1 )), 



a(G 2 ,H 2 )=(3(a(G 1 ,H 1 )). 
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In this section and the next one, we show that o~(G, G) is a large subset of 
A(G,G), in the sense that it contains a linear subspace of co-dimension < 5. Con- 
sequently, by the existence of bijection (3, the same must be true for a(G, H) inside 
A(G, H) when G and H are locally equivalent. This observation suggests us to 
consider just one graph instead of two, i.e. to assume that G = H. 

Definition 4.1 Internal solutions for a graph G are vectors (X, Y, Z, T) in a(G, G). 
4.1 A(G, G) 

The orthogonality assumption, equation ([2]), in the case of G = H can be written 
more efficiently. Indeed, it is easy to check that (|2|) is equivalent to 

(X,g(i) x g(j)) = (Y(j) - Z(i)) 9ij , for every i + j, (3) 

(X,g(i) x g(i)) = T(i), for every i. (4) 

Lemma 4.1 Assume that the graph G is connected. Then, for every (X,Y, Z,T) 
£ A(Gr, G), the function Y + Z on the vertices is constant, i.e., Y + Z = al for some 
number a. 

Proof: Assume that i,j are two adjacent vertices in G. The condition ([3]) implies 
that 

(X,g(i)xg(j)) = (Y(j)-Z(i))g ij7 

(X,g(i)xg(j)) = (Y(i)-Z(j)) gji . 

Therefore Y(i) + Z(i) = Y(j) + Z(j), for any two adjacent vertices. Connectivity 
of G implies the desired conclusion. 

□ 

From now on, we assume that G is a connected graph. Lemma 14.11 gives us a 
partitioning of the set A(G, G), as follows: 

A a (G, G) = {(X, Y, Z, T) G A(G, G) :Y + Z = al}. 

Notice that, for any a E F q , 

(0, al, al, 0) G A 2a (G, G), 

and 

A 2a (G, G) = Ao(G, G) + (0, al, al, 0). 
Therefore, since q is an odd number, all A a (G, G)'s are just shifts of Aq(G, G). 
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Definition 4.2 Consider a connected graph, G. We say that ij is an edge of G if 
gij ^ 0. For an even cycle C in G consisting of (ordered) vertices i\, 12, ■ ■ ■ 121, 1st 

21 

< C ) = D" 1 ) k 9r k \ k+ Mik) x g(i h+1 ). 

k=l 

We define v{G), called the bineighborhood space of G, to be the subspace generated 
by the vectors g(i) x g(j) for i,j satisfying g^ = 0, as well as by v(C) 's for even 
cycles C, i.e., 

v(G) = Span{{v(C) : C even cycle) U {g(i) x g(j) : g^ = 0}}. 



Theorem 4.1 For X G F£, there exist Y, Z and T such that (X, Y, Z, T) G A (G, G) 
if and only if X is orthogonal to v{G). Moreover, if G has an odd cycle, for every 
X G v{G) L , there exists a unique (X,Y, Z,T) in A (G,G). 

Proof: Fix a vector X, and assume that there exists some (X, Y, Z, T) G Ao(G, G). 
For any two vertices i,j with = 0, (X,g(i) x g(j)} = (Y(i) +Y(j))g i:j = 0. 
Moreover, if ii,i2, ■ ■ ■ ^21 is a cycle then we have 

(X, (-lfg-^ +i g{i k ) x g(i k+1 )) = (-l) k (Y(i k ) + Y{i k+l )). 

By summing up all of these equalities for k = 1, 2, • • • , 21 we get that X is orthogonal 
to f(C), for every even cycle C, and hence, to the whole space v{G). 

For the other direction, suppose that X G v{G)- i ~. For any vertex i, set T(i) = 
(X,g(i) x g(i)), and Z = —Y. In the case that G has an odd cycle, Y would be 
uniquely determined by ([3]) on the vertices of that odd cycle. Since, G is connected, 
then one can determine the function Y on the rest of the vertices, and since X G 
v(G)~ L there is no ambiguity in the definition of Y. 

In the case that G contains no odd cycles, we can fix Y(i) for some arbitrarily 
chosen vertex i, and then determine the other components in terms of Y(i). Once 
again, since X G ^(G) 1 - and there is no odd cycle, there is no ambiguity in its 
definition. 

□ 

4.2 Vectors in A(G, G) have constant determinant 

Even though, the determinant is a quadratic function and not a linear one, but in 
this setting, the set A(G, G) satisfies some property that helps us to study this set 
more deeply. 

The problem of verifying locally equivalent graphs on one or two vertices is a 
trivial problem, and hence, from now on we assume that the number of vertices of 
G is more than 2. 

Theorem 4.2 For every (j> G A(G,G), the determinant of <f> is constant. 
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Proof: First, notice that one may restrict oneself to the case <ft G Aq(G,G), 
since there exists a vector (X,Y, —Y,T) G Ao(G, G) and a G F g such that (ft = 
(X,Y,-Y,T) + a(0,I,I,Q), and det<ft = det(X,Y,-Y,T) + a 2 I. Hence, suppose that 
4> = (X, Y, — Y, T) G A (G, G). 

If G contains at most two vertices, the proof is clear. So, assume that n > 3, 
and ii,i 2 are two adjacent vertices in G. Showing that (detcft)^ = (det<p)i 2 gives 
us the desired result. By lemma |3T2| local complementing operations, * and o, do 
not change the determinant, therefore we can assume that i\ and i 2 have a common 
neighbor j$ by applying one * operator if needed. One has 

g rs (Y(r)+Y(s)) = (X,g(r) x g(s)), 

for each pair of unequal r, s G {jo, ii, 22 }• Consequently, 



l r 



= 2[4( 1 '^) x sfe)) + <7(»l) x g(j Q ))-gr) o {X,g(i 2 ) x 5 (j )) 



and 



^2) = ^[gr i U X ^(k) x 9 (i 2 )> + 54<X, 5 (i 2 ) x g(j ))-g£ (X,g(h) x 5 (j )> 



On the other hand, T(r) = (X,g(r) x g(r)), for any vertex r. Hence, in order to 
prove (detcft)^ = (det(ft)i 2 , we should show that 



-l[9hl( X ,9(h) x g{i 2 )) +gr) o {X,g(i x ) x s (j )) 
-g-l(X,g(i 2 ) x g(j ))] 2 - X^X^ih) x g(h)) 
= -4 ^i 2 ( x >5(n) x 5 (« 2 )} +^) o (X, 5 (i 2 ) x 5 (j )) 
-^(X^i) x 5 (j ))] 2 - X(i 2 )(X,g(i 2 ) x 5 (i 2 )>, 

or equivalently 

^! 2 (^5(n) x g(i 2 )) g^X^fa) x g(j )) - g^{X, g(h) x 5 (j )) 
= Xii^gih) x g^)) - X(i 2 )(X,g(i 2 ) x g(i 2 )). 



Let 



Cj = g-){X,g( i2 ) x g(j)) - g^(X,g(h) x g(j)), 



for any j adjacent to both %\ and i 2 . Since, X is orthogonal to the cycle jo,ii,j,i 2 , 
for any j adjacent to ii,i 2 , we have Cj = Cj. On the other hand, if either g^j or 
gi 2 j is zero, then 

9iij9i2j x {j)Cj = 0, 
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and 



X{3){9hj{ X ,9(i2) x 9{j)) -9i2j{ x ,9{h) X g{j))) = 0, 



because, for instance if g^j = 0, then X is orthogonal to g(i\) X Therefore, we 
have 



(X,g(ii) x 9{h))[g i2 ) (X,g{i 2 ) x 9(h)) -gi^ {X,g(h) x 5 (j )) 
Y 9hj9i 2 j x U)) 9i 2 j (X,g(i2) x <?(j )} - g^iX^ih) x g(j )} 

Y 9i 1 j9i 2 j x U) 9i 2 j ( x ,9(i2) x g(jo)) ~ g-^iX^ih) x g(j )) 

Y 9i 1 j9i 2 j x U) [9rjj( x >9(i2) x g(j)) - g^iX^ih) x g(j)) 

Y X W ghi( x ,g{n) x g(j)) - gi 2 j{x,g{h) x g (j)) 

ii 

= Y Y 9i d9i 2 k9kj x (j) x (k) - 9i 2 jgi 1 kgkj x (j) x (k) 

= + Y Y 9hjgi2kgkjX(j)X(k) - g i2j g ilk g kj X(j)X(k) 
j^ii,i 2 k=h,i 2 
n n 

= Y 9 hi2 9i 1 j X ( i i)X(j) - Y9hi2gf 2 j X ( i 2) x U) 

g ili2 ^X(ix){X,g(ii) x g(i x )) - X(i 2 )(X,g(i 2 ) x g(i 2 )) , 
which completes the proof. 

□ 



5 Linearity of the kernel of det function 

Using theorem 14.21 we may give another partition of the set A(G, G) as follows, 

A a (G, G) = {<p£ A(G, G) :det<\> = al}, 

and combining with the previous partition, we set 

K{G,G) = A a (G,G)nA a (G,G). 

Notice that Ao is a linear subspace since it is the kernel of a linear map. But 
generally the determinant function is not a linear function when q is not a power of 2. 
Here, due to the nature and strength of theorem !4.2l we will show that despite of not 
being a linear function, the kernel of the determinant exhibit some linear properties. 
More precisely, we show that A[j(G, G) is a linear subspace if dimA(G,G) > 5. 
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5.1 Some useful lemmas 

For any 

= (X, Y, -Y, T), 4> = (X', Y', -Y', T') G A (G, G), 

define 

0') = 2Y x Y' + X x T' + X' x T. 
Notice that, ^(0, <j)') is constant, because det is a constant function and 

0') = det + det (p f - det (0 + 0'). 

Also, to prove the linearity of A[](G, G), one can sufficiently show that Vl/(0, fa) = 
for any (p,fa G Aj] (G,G). The following series of lemmas provide us the necessary 
tools. 

Lemma 5.1 Suppose that fa = (Xj, Y, — Y, Tj) G Ao(G, G) /or i = 1,2, and more- 
over, fa G A[|(G, G). suppose that ^(fa,^) = al for some 0^o£ F q . Then 
on every vertex, either X\ or X<i is non-zero. 

Proof: Assume that X\(i) = 0, for some vertex i G {1,2, • • • , n}. Since det (pi = 
0, one has —Y\(i) 2 — Xi(i)T\(i) = and therefore, Y\(i) = 0. Hence, the i-th 
component of ^(fa, (j>2) is X2(i)T\{i) = a/ and thus X2(i) ^ 0. 

□ 

Lemma 5.2 Let fa = (X h Yi,-Yi,Ti) G t§{G,G) fori = 1,2, andiP = (U,V,-V,W) 
G A[](G, G). Moreover, assume that ^>(fa, $2) = al for some / o £ F q , and 
^(4>i, ip) = 0. Then supp(U) C supp(Xi), in the sense that if U is non-zero on 
some vertex, then so is X\. 

Proof: For any r G F q , we have \P(0i,r^> + $2) = al ^ 0. Therefore, by 
lemma [oTTI supp(X\) U supp(rll + X2) = {1, 2, • • • , n}. Now suppose that X\(i) = 
for some % G {1,2,- ■■ , n}, and f7(i) 7^ 0. There exists some ro G F g such that 
roU(i) + X2(i) = 0. Thus « ^ supp(Xi) U supp(roU + X2), which contradicts the 
earlier statement. 

□ 

The third lemma is the following: 

Lemma 5.3 Suppose that fa = (Xi,Y{, —Yi,Ti) G A[](G, G) for i = 1,2, such that 
supp(X\) and supp(X2) are minimal subsets o/{l,2, ••• , n} with ^(fa,^) 7^ 0. 
IfiP = (U,V,-V,W) G Ag(G,G) and = 0, i/ien e^/ier ^ = or [/ = 

— a _1 cXi, where ^{fa,4>2) = o-f ana 1 ^(^VO = c -^- 

Proof: First, assume that c = 0. In this case we have ^(rip + fa,4>2) = o,I, for 
any r G F 9 . If J7(z) 7^ for some i G {1, 2, • • • n}, then there exists ro G F q such 
that roU(i) + X\(i) = 0. Thus, using lemma [572], supp(U) C s-upp(Xi) and hence 
supp(roU + Xi) is a proper subset of supp(Xi). Moreover, ^(roV* + fa, $2) = o,I 
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and det(rQtp + fa) = rfidet ij> + det fa - r ^(ip, fa) = 0. Then r ip + fae A° Q (G, G), 
which contradicts the minimality of fa, fa. Therefore U = 0. 

Now, assume that c / 0. Once again, by lemma supp{U) C supp(Xi) 
and supp(Xi) C supp(U). Suppose that U(i) ^ for some i G {1, • • • , ra}. There 
exists some r G F g such that r U(i) + X^i) = 0, and also, ^(rQip + ^1,^2) = 
(roc+a)I. By the minimality of fa and 02) one concludes that r^c+a = 0. Therefore, 
—ac~ l U{i) + Xi(i) = for any i satisfying U(i) 7^ 0. Therefore U = —a~ l cXi. 

□ 

Finally, the last lemma is a fact in number theory. 

Lemma 5.4 Given a 3 x 3 matix A in F q , there exists a non-zero vector x g so 
that x T • A ■ x = 0. 

Proof: There are many ways to prove this lemma, and one straight forward com- 
putational way is as follows. Let us rewrite the matrix equation as a degree two 
numeric equation, 

ax 2 + by 2 + cz 2 + 2dxy + 2exz + 2fyz = 0. 

If abc = then there exists a trivial solution to the equation. For instance, (1,0, 0) 
when a = 0. Thus, we assume that abc 7^ 0. 

Solving the latter equation in terms of z using the square root of delta formula, 
we obtain that the problem is equivalent to this one: does delta have a square root? 
In other words, it is equivalent to finding a non-trivial solution to the following 
equation, 

ax 2 + [3y 2 + 2jxy = t 2 , 

where a = e 2 — ac, j3 = f 2 — be and 7 = ef — cd. Once again, if a = then 
(x,y,t) = (1,0,0) is a solution, and if not, by solving it in terms of x, we get the 
next equation, 

y 2 = s 2 - at 2 , 

where 9 = -y 2 — a/3. We set y = 1. For different values of s,t G F q , each of the 
functions s 2 — 8 and at 2 ranges over (q + l)/2 different elements. Therefore, there 
is at least one (s,t), such that s 2 — 9 = at 2 , and we are done. 

□ 



5.2 A° (G, G) is linear 

Now, we have all the necessary tools in hand, to provide a proof for the linearity of 

Theorem 5.1 Suppose that dimA (G, G) > 5. Then ^ = on A(j(G, G), or equiv- 
alently, A[j(G, G) is a linear subspace. 
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Proof: First of all, we can assume that G has an odd cycle. Because we already 
know that by local complementation, the linear properties of A(G, G), the determi- 
nant and so the function do not change. Under this assumption, as we observed 
in theorem SH if <fc = (Xi, Y i: —Y{, Tj) € A (G,G), i = 1,2, and Xi = X 2 , then 

01 = 4>2- 

Suppose that, ^ is not zero, and let (pi = (Xi, Yi, —Yi,Ti) € A[](G, G), i = 1,2, 
such that Xi, X2 are minimal elements of {1, 2, • ■ ■ , re} (in the sense of lemma [5T3]) 
satisfying ^(0i,02) = al, where ^ a € F g . Since, dimA^G, G) > 5, there 
exist elements -0? = (Uj,Vj, —Vj,Wj) € A (G, G), j = 1,2,3, independent of 0i 
and ^2- Set ^((pi,tpj) = bj and ^((p2-,ipj) = Cj for j = 1,2,3, and also define 
iv j = aipj — Cj(pi — bj4>2, for j = 1, 2, 3. One can easily verify that ^f((pi, u>j) = 0, for 
every i = 1,2 and j = 1,2,3. Using lemma [£T4"1 we can find a non-trivial solution 
of det (r\U)\ + r 2 W2 + r 3 uj 3 ) = 0, where n, r 2 , r 3 G F g . Thus, ri^i + r 2 ^2 + ^3^3 € 
Aj](G, G) and ^(^jriwi + r2UJ2 + ^3^3) = 0, i = 1,2. Therefore, by lemma f5T3l the 
first coordinate of r\LO\ + ri^2 + r 3^3 is zero and also by theorem 14. \\ we conclude 
that r\<jj\ + T2U)2 + ^3^3 = 0, which is a contradiction. Hence, ^f(4>i, $2) = for 
every (pi,(p2 £ A[j(G, G), and A[j(G, G) is a linear subspace. 

□ 

This theorem says that Aq (G,G) is a linear subspace. But, we can say much 
more about that. Having constant determinant as well as its linearity, makes it 
a significantly helpful to study o~(G, G) whose description is our main goal in this 
section. In fact, the following lemmas tell us that this linear space, A[j (G, G), is 
really a large subspace in the whole space A(G, G). 

Lemma 5.5 The co-dimension o/A[](G, G) in Aq(G,G) is at most two, provided 
that dimA (G, G) > 5 . 

Proof: Consider three independent vectors (pi, (p2,<Pz £ Ao (G,G). For numbers 
ci,c 2 and c 3 , 

det (c\(pi + c 2 (p 2 + C303) = c\det ((pi) + c\det (0 2 ) + c\det ((p 3 ) 

+CiC2^((pl,(p 2 ) + Cl C 3 #((pl, (p 3 ) + C 2 C 3 ^((p 2 , 03). 

By lemma [5^1 there exists (ci, C2, c 3 ) 7^ such that det (ci0i + C20 2 + £3^3) = 0, 
which means that ci(pi+C2(p2+c 3 (p 3 G A{j(G, G). Thus, the co-dimension of A[j(G, G) 
inside Ao (G,G) is at most two. 

□ 

Since q is an odd number, the translation of Ao(G, G) by vectors (0,al, al, 0), 
for different values a G F g , gives the whole space A(G, G). Therefore, co-dimension 
of A (G,G) in A(G,G) is one. Also, Ag(G,G) + (0,1,1,0) = A|(G,G). On the 
other hand, by the definition of cr, <r(G, G) 5 A^G, G). Therefore, we conclude the 
following corollary: 

Corollary 5.1 There exists an affine linear subspace inside o~(G,G), whose co- 
dimension in the whole space A(G, G) is at most 3, provided that A(G, G) has di- 
mension not less than 5. Putting these together, in general the co-dimension of this 
affine subspace is at most 5. 
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Having the mentioned property in hand, together with the (3 function introduced 
earlier, we can now give the desired description of o~(G, H ) for equivalent graphs G 
and H, which creates the foundations of our algorithm determining whether two 
graphs are equivalent or not. 

Theorem 5.2 If the connected graphs G and H , defined on the same vertex sets, 
are locally equivalent, then the co-dimension of some affine linear subset of a(G,H) 
inside A(G,H) is at most 5. 

6 The algorithm 

We now have all the tools to describe in details and provide the proof for an efficient 
algorithm to determine whether two graphs are equivalent or not. The algorithm is 
the following. 

Suppose that G and H are two connected graphs (notice that by local comple- 
mentation a connected graph remains connected), with neighborhood functions g 
and h. Consider the linear system of equations: 

(X,g(i) x h(j)) - (Y,g(i) x e 3 ) + (Z, ei x h(j)) - (T, e % x e 3 ) = 0, (5) 
for any two vertices and for X,Y, Z and T in F™, together with the equation 

Y x Z - X x T = I. (6) 

Assume that B is an arbitrary basis for A(G,H), the set of solutions of the linear 
equation ([5]), which can be computed efficiently. According to the corollary 1 5 . 1 1 and 
theorem I5.2l in the previous section, if there exist solutions for (JSJ) and ([6]), then there 
exists an affine subset, denoted by T, in A(G, H) with codim < 5, whose elements all 
satisfy both ([5]) and (|6|). The following lemma takes the advantage of this property. 

Lemma 6.1 For any basis B of a linear space A, and every affine subspace T of A 
of codim < 5, there exists a vector u G T, which is a linear combination of at most 
five elements of B. 

Proof: Consider the set V = {u — v : u, v S T} which is a subspace of A, 
and the canonical projection p : A — > A/V . The set p(B) generates A/T', and 
therefore there exists a basis {p(b\), . . . ,p(bk)} for A/T' , where b\,...bk £ B and 
k = dim (A) — dim (V) < 5. Since T is affine, T G A/T' and hence can be written 
as the linear combination of p(6j)'s, which means that there exists a vector ueT 
which is a linear combination of at most five elements of B. 

□ 

By this lemma, we can now consider all of the linear combinations of every 5 
elements of B, and check whether or not it satisfies the condition ([6]). If at least one 
of them satisfies ([6]), then the answer is positive, and is negative otherwise. 

Notice that, solving this problem for disconnected graphs is an immediate con- 
sequence of solving it for the connected graphs, since the local operators preserve 
the connectivity. 
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The described algorithm is efficient. In fact, by using a pivoting method, a basis 
B can be computed in 0(n 4 ) time, because there are 0(n 2 ) linear equations in ©. 
This number must be added to and hence will be dominated by the time to check 
the equation ([6]) for all of the linear combination of five elements of this basis, which 
is 0(n 5 ), (in the case that dimA(G,H) < 5, we check all of the possibilities). Thus, 
the algorithm takes 0(n 5 ) time, and the overall complexity is polynomial in n. 
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